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ABSTRACT 

Using adaptive optics assisted Gemini/NIFS data, I study the present and past gas accretion in 
the central 3" of the M32 nucleus. From changes in the spectral slope and CO line depths near the 
center, I find evidence for unresolved dust emission resulting from BH accretion. With a luminosity 
of ~ 2 x 10 38 ergs -1 , this dust emission appears to be the most luminous tracer of current BH 
accretion, two orders of magnitude more luminous than previously detected X-ray emission. These 
observations suggest that using high resolution infrared data to search for dust emission may be an 
effective way to detect other nearby, low luminosity BHs, such as those in globular clusters. I also 
examine the fossil evidence of gas accretion contained in the kinematics of the stars in the nucleus. 
The higher-order moments (/i3 and M) of the line-of-sight velocity distribution show patterns that 
are remarkably similar to those seen on larger scales in elliptical galaxies and in gas-rich merger 
simulations. The kinematics suggests the presence of two components in the M32 nucleus, a dominant 
disk overlying a pressure supported component. I discuss possible formation scenarios for the M32 
nucleus in the context of the kinematic data as well as previous stellar population studies. The 
kinematic measurements presented here are the highest quality available for the nucleus of M32, and 
may be useful for any future dynamical models of this benchmark system. 

Subject headings: galaxies:nuclei - galaxies: active - galaxies: kinematics and dynamics - galax- 
iesdormation - galaxies: individual (M32) 



1. INTRODUCTION 

The M32 nucleus hosts the smallest of the three super- 
massive black holes known in the Local Group. Like its 
counterpart in M31 and the Milky Way, it is currently 
accreting at very low levels, detectable only due to its 
proximity. Surrounding the black hole (BH) in M32 is a 
rotating disk of stars that is the densest stellar system 
in the local universe. In this paper I use high resolution 
spectroscopic observations to examine the current accre- 
tion of the M32 black hole and study the formation of 
the stellar nucleus from its present-day kinematics. 

The kinematic signature of a super massiv e BH a t the 
center of M32 was first recognized by Tonry (1984), and 
current estimates su ggest a BH mass of 2.5xl 6 M Q 
(jVerolme et al.1 12002t1van den Bosch fc de Zeeuwll2010h . 
Accretion onto this BH has been searched for at many 
wavelengths, but detected only in the X-rays (|Ho et al.l 
I2003h . The X-ray luminosity from 2-10 keV is 9.4xl0 35 
ergs -1 , suggesting an accretion rate nine orders of mag- 
nitude below the Eddington limit. Like M31 and the 
Milky Way, the M32 BH presents a rare opportunity to 
study BH accretion processes at very low levels. In §3, I 
report detection of hot dust emission from near the black 
hole at NIR wavelengths with luminosity exceeding the 
X-ray luminosity by two orders of magnitude. 

In the second half of the paper (§4), I use the kine- 
matics of the nuclear region to constrain the formation 
history of its stars. The formation history of stars in 
the immediate vicinity of the BH is likely to be con- 
nect ed in some way to the form ation of the BH itself 
(e.g. iHopkins fc Quataerd l2010bl) . The nucleus of M32 
contains a rotating disk of stars that reaches a den- 
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sity at r < 0.1 pc of >10 7 M^/pc" 3 , making it the 
dense st stellar system known (|Walkerlll962l : lLauer et al.l 
,1998:). Morphologically, the nuclear component appears 
to be distinct from the rest of the galaxy, with the 
surface brightness profile showing evidence f or an addi- 
tiona l stellar component at rad ii <5" (20 pc) (|Kormendvl 
119991: IGraham fc Spitlerl 120091 ) . This component has an 
effective radiu s of 1.5", and an I b and luminosity of 
-3xl0 7 M (IGraham fc Spitlerl 12001 . ~10% of the to- 
tal galaxy luminosity (see also iKormendv et all [2009). 
The mass and size typical of this component is typical 
of nuclear star clusters commonly seen at the centers of 
lower luminosity early-type galaxies, but it contains a 
much larg er fraction of the total galaxy luminosity than 
usual (e.g. lCote et al.ll200l . 

Despite the evidence for a morphologically distinct 
component in the central 5", no corresponding sepa- 
ration is seen in the stellar populations or kinemat- 
ics. The rotation velocity d ecreases slowly with radius 
([Dressier fc Richstond Il988) but there is no break sug- 
gesting a separation between the nuclear disk and the rest 
of the galaxy. The stellar populations of the nucleus are 
younger, than the galaxy as a whole, with an average age 
of ~4 Gyr in the nucleus rising to ~8 Gyr at larger rad ii 
(|Worthevl [2001 IRose et al.l l200l ICoelho et all 12009ft . 
However, again there is no evidence for a break in the 
stellar population properties as a function of radius. In 
§4 I show that there is a kinematic signature in the nu- 
cleus that is remarkably similar to merger simulations 
done at larger scales. This signature provides strong ev- 
idence for a two component system in the M32 nucleus, 
with a dominant disk component embedded in a pressure- 
supported component. I outline a scenario, in which the 
nuclear disk formed gradually from the stellar winds of 
stars in the bulge of M32, that may provide a good ex- 
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planation of the observed kinematics, stellar populations 
and abundance gradients seen in the nucleus. 

2. DATA 

Integral field spectra of the central 3" x 3" of M32 were 
obtained with Gemini NIFS on Oct 23, 2005. A total of 
11 x 600s K band exposures were taken, with 6 on source 
and 5 offset to a blank sky position. The nucleus was 
used as a natural guide star to obtain an adaptive optics 
correction with the ALTAIR syste m. Two previous pa- 
pers have utilized this same dataset (|Davidge et al.ll2008l 
2010') but were restricted to the study of stellar popula- 
tions in the nucleus. 

Data were obtained from the Gemini Science Archive 
and reduc e d fol lowin g the sam e proc ess outlined in 
ISeth et alJ (12008ft and ISeth et all ([2010ft . An A0V star, 
HIP 116449, was used as a telluric calibrator. This 
star was also used for photometric and spectrophoto- 
metric calibration, both of which have errors of ~10%. 
Wavelength calibration utilized both an arc lamp image 
and the sky lines in the spectrum, and has an abso- 
lute error of ~2 kms -1 . The median spectral resolu- 
tion was 4.2A (57 kins" 1 ) FWHM. The signal-to -noise 
ratio of the spectrum is very high, ranging from 280 per 
0'.'05" x 0'.'05 pixel in the center to 60 in the corner of the 
field (r = 2'.'l). 

To determine the point spread function (PSF) of the 
observations I followed a two-st ep process similar to that 
described in ISeth et all (|2010l ). I assumed a PSF de- 
scribed by an inner Gaussian core + outer Moffat halo 
(with the Moffat function described as E(r) = So/[(l + 
{r/rd) 2 } 4 ' 765 )- I used the observation of the telluric cali- 
brator to measure the outer halo of the PSF and found 
a best fit = 0'.'83 containing 55% of the light. I then 
fixed this outer halo, and fit the inner part of the galaxy 
PSF by convolving an HST image of M32 to fit our NIFS 
image (after making a minor correction for the dust emis- 
sion, see §3 for details). For the HST image, I used an 
F1042M WFPC2 image where the nucleus was located 
on the planetary camera. I found a very good fit to the 
NIFS image by convolving this HST image using an in- 
ner Gaussian with FWHM of 0'.'25 and containing 45% of 
the light (along with the outer Moffat with rd — 0'.'83); 
typical fit residuals were 4%. The compact dust emission 
discussed in the next section has a FWHM of 0'.'25, an 
independent confirmation of the PSF core width. Mea- 
surements of the CO bandhead strength used in §3 and 
plotted in Fig. [T] were made in each spaxel after dividing 
by the tell uric star spectrum using t he 12 CO (2, 0) index 
defined bv lKleinmann fc Halll (|1986ft . 

I derived kinematics from the NIFS data from 
the CO bandhead between 2.28/^m and 2.395/xm us- 
ing the penalized-pixel fitt ing (PPXF) method of 
ICappellari fe Emselleml (|2004ft to determine the line-of- 
sight velocity distribution (LOSVD). The program pa- 
rameterizes the LOSVD using the radial velocity (V), 
dispersion (tr), skewness (h3) and kurtosis (fe4). High 
resolution templates from lWallace fc Hinkld (|1996ft were 
used and convolved to the spectral r esolution of each 
pixel as determined from sky lines (see ISeth et al.ll2"010L 
for more details). 

The measurements presented here are the highest qual- 
ity kinematic data available that resolves the sphere- 
of-influence of the BH. The very high S/N results in 
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Fig. 1. — The CO line strength plotted against continuum slope. 
Positive values indicate deeper CO lines and redder continuum 
slopes. Spaxels between O'.'l and 0'.'2 of the center are shown with 
blue crosses, and those within O'.'l with red crosses. The deviation 
towards redder colors and weaker CO lines is due to emission from 
hot dust. The slope is measured from the flux ratioed spectra, thus 
a value of zero indicates the same spectral slope as an AO star. Er- 
rors on the CO depths are <0.005 and on the continuum slopes, 
<0.0002. 

errors on the LOSVD that are much sm aller than for 
the H ST/STIS measurements presented bv lJoseph et alJ 
(|2001[) . The spatial resolution of the NIFS observa- 
tions is somewhat lower than with STIS, but NIFS pro- 
vides integral field coverage that can hel p significantly 
in constraining dynamical models (e.g. [Shapiro et alJ 
12001 iNeumaver e'taLI 120071 iNowak et all 12008ft . Our 
kinematic m e asure ments are consistent with those of 
IJoseph et al.1 (|2001|) except in the central few tenths of 
an arcsecond, where the higher resolution of STIS re- 
sults in a steeper velocity gradient and a central dis- 
persion higher by ^10%. As M32 is a benchmark 
for stu dies of BH mass measurements and dynamics 
codes (Ivan der Marel et all 119981: iVerolme et all 120021 : 
iKormendvl 12004 Ivan den Bosch fc de Zeeuwl I2010Q . I 
make the full kinematic data set available in Table 1 for 
use in future studies. 

3. HOT DUST: THE MOST LUMINOUS TRACER OF 
CURRENT BH ACCRETION 

The spectral energy distributions (SEDs) of quasar 
and Seyfert nuclei in the K band are typically dom- 
inate d by hot dust with temperatures of 800-1500 K 
(e.g. H Cobayas hi et al.l 119931: lAlonso-Herrero et al.l 119961 : 
IWinge et all 120001: IRiffel et all 12009ft . High resolution 
observations with HST NICMOS and adaptive optics ob- 
servations have shown that these sources are typically un- 
resol yed, and thus originate quite close to the accreting 
BH (|Quillen et al.|[200lt IRiffel et al.l [2010h . From mea- 
surements of the time-lag between UV and NIR light, 
the size of this dust-emitting region has b een constrained 
to be -Cljpc in nearby Seyfert galaxies (|Minezaki et al.1 
12001 12001 . The luminosity of the dust emission cor- 
relates s trongly with AGN luminosities at other wave- 
lengths ([Guillen et al.1 12001ft . Recently, we found unre- 
solved and possibly variable hot dust emis sion in the 
nucle us of the nearby SO galaxy NGC 404 (jSeth et al.1 
12010ft . the first detection of compact dust emission in a 
low-luminosity (LINER) AGN. 

The presence of dust emission in the nucleus of M32 
is inferred due to the presence of a redder continuum 
and weaker CO lines in the central 0'.'2. This is shown 
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TABLE 1 

Gemini/NIFS Kinematic Data of the M32 Nucleus 



Aa ["] A8 ["] V r [kms" 1 ] <r [kms" 1 ] h3 hA 

0.014 -0.015 -206.14±1.08 117.94±1.23 0.018±0.007 0.038±0.007 
0.014 0.035 -215.93±1.09 114.63±1.17 0.022±0.007 0.053±0.007 
-0.036 -0.015 -208.52±1.10 115.40±0.87 0.022±0.005 0.043±0.007 
(See Online Journal Article for Full Data Table) 



in Fig. [TJ Typically, the color of the continuum corre- 
lates with the depth of the CO line as cooler redder stars 
have deeper CO lines. However, within the central 0'.'2, 
(and especially the central O'.'l shown with red crosses 
in Fig. [TJ), the continuum gets redder while the CO lines 
become weaker indicating the presence of hot dust emis- 
sion. To determine the contribution from this hot dust 
component I fit the M32 nuclear spectrum to a combina- 
tion of a blackbody model an d an annular spec trum from 
radii of O'.'5-l" (as in Fig. 4 of lSeth et al.l[2010l) . The best 
fitting temperature for the dust at the center (r < O'.'l) 
is 920K. The residuals of the fit are 1.1%, a factor of two 
improvement over a fit with no dust. 

I fit the dust contribution to each spaxel in the data 
cube fixing the dust temperature to 920K. Dust emission 
is only found at r <0'.'2 and reaches a maximum contri- 
bution of ^5% of the total light. The dust emission is 
unresolved, with a FWHM of 0'.'25 matching the NIFS 
PSF derived from HST observation. The emitting region 
is therefore constrained to be <0.9 pc. 

The total flux of the hot dust emission from 2.0- 
2.43 /mi is 1.7xl0 -13 ergs -1 cm -2 corresponding to a 
luminosity in the K band of 1.2xl0 37 ergs -1 . Integrat- 
ing over the full Planck function, the total luminosity 
is 1.8xl0 38 ergs -1 . The corresponding K band mag- 
nitude of the emission is if=15.0 (M K = -9.4). The 
primary uncertainty in these flux and luminosity mea- 
surements is the temperature of the dust, which is not 
well constrained by our fits. Taking a reasonable range 
of temperatures from 800-1300K, the flux ranges from 
1.5-2.5 xlO -13 ergs -1 cm -2 , the K band luminosity from 
1.0-1.7xl0 37 ergs -1 and the total luminosity from 1.3- 
2.6xl0 38 ergs -1 . 

The weaker CO lines and redder continuum cannot re- 
sult just from changes in the stellar population at the 
center of M32. There is little indication of any change 
in the optical and UV colors near the center of M32, al- 
though slightly bluer colors are seen in the c entral pixels 
of a W FPC2 F555W-F814W color image bv lLauer et all 
(1998). The Galactic center also sh ows a decrease in CO 
line depth within the central parsec ( Scllgrcn "et~aT1[T990l) 
which likely results from the decrease in the fraction of 
light coming from l ate-type stars i n the central parse c 
(|Genzel et al.llT99l IPolrtaTl l200l iBartko et al.ll2010h . 
This decrease may be caused by stellar collisions pref- 
erentially destroying the la rger radius RGB and AGB 
stars (e.g. iDale et all |2009f ). The M32 stellar density 
is sufficient that s uch collisions may occur there as well 
(jLauer et alj|l99 8). However, this mechanism cannot by 
itself explain the effect shown in Fig. (TJ as the removal of 
larger radius RGB and AGB stars would make K band 
continuum bluer, not redder. 

Emission from hot dust could also result from heating 
by the stellar component. M32 has a population of post- 



AGB and ex treme Horizontal Branch stars ([Brown et al.l 
[19981 l2000l) . but these are much less compactly dis- 
tributed than the dust emission. Furthermore, the total 
UV luminosity at 1600A is insufficient to heat the dust 
{vL v - 2 x 10 37 ergs -1 ; ICole et al.l[l998h . and the UV- 
to-optical colors d o not decrease at the center of M32 
(jLauer et al.|[l99 8j). Therefore hot stars do not appear 
to be responsible for the dust emission. The interact- 
ing stellar winds resulting from the high stellar density 
and velocity dispersion could also provide energy to heat 
the dust. The conditions at the center of M32 are fairly 
similar to those in the MW center where no comparable 
hot dust component is seen. However, there is a signifi- 
cant warm dust component in the Milky Way center with 
T~400 K and lu minosity of order 10 38 ergs -1 within 
the central parsec (|Zvlka et al.lll995HStolovv et al.lll99l 
IMezger et al.lll996l ). I also note that there are no signa- 
tures of compact hot dust emission in NIFS observations 
of local group galaxies M33 and NGC 205 which also have 
high stellar densities (but lower dispersions). Based on 
the lack of hot dust emission in these systems, I argue 
that the hot dust in the nucleus of M32 is likely powered 
by BH accretion. 

The compact hot dust emission appears to be the most 
luminou s BH accretion signature observed in the M32 
nucleus. IHo et al.l f|2003T ) analyzed a variety of radio and 
X-ray data and detect the nucleus only in X-rays, with a 
2-10 keV luminosity of 9.4xl0 35 ergs -1 . Based on this 
X-ray lumin osity and the exp ected supply of gas available 
for accretion IHo et ahl ([2003) found that the radiative or 
BH accretion efficiency must be very small (< 10 -3 ). 
Such low efficiencies a re found in radiatively ine fficient 
accretion models (e.g. iQuataertl |200H iHol 120091 ). The 
NIR emission discussed here would allow for a somewhat 
higher accretion efficiency, but still likely lower than the 
10% efficiency of a canonical accretion disk. This sug- 
gests that strong NIR dust emission may be found in 
lower accretion states as well as the higher accretion 
states where they have traditionally been observed. The 
dust emission in Seyferts and quasars is thought to come 
from a radiatively heated sublimation zone in the accre- 
tion disk or at the inner edge of the torus. However, the 
excitation mechanism for the hot dust emission in M32 is 
unclear. The NIR luminosity derived here is two orders 
of magnitude brighter than th e observed X-ray luminos- 
ity found by IHo et all (|2003l ) . In studies of other low 
luminosity AGN , the NIR to X-r ay luminosity rat ios are 
typically unity (|Hol 119991 120081 : lEracleous etall l2010f ). 
These systems are somewhat higher luminosity systems 
than M32, but all are expected to be radiatively ineffi- 
cient accr etors with similar spectral energy distributions 
(|Hol [2009) . Furthermo re, the observed X -ray luminosity 
and UV upper limits (|Cole et al.l Il998t IHo et al.l 120031) 
are insufficient to supply the radiative energy needed to 
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Fig. 2. — The line-of-sight velocity distribution of M32 derived from the Gemini/NIFS data at the CO bandhead (2.3/^m). The radial 
velocity is shown in the top left panel after removal of the systemic velocity of -200 kms -1 . The bottom left panel shows the dispersion, 
the top right the skewness or h3 component, and the bottom right the kurtosis or h4 component of the LOSVD. These data are the highest 
S/N data available of the M32 nucleus that resolve the BH sphere of influence. 



heat the dust. Both these problems could be solved if the 
BH accretion in M32 is significantly time variable, with 
our observations taken at a time when there was higher 
flux at all wavelengths. Alternatively, the heating of the 
hot dust could result from m echanical feedback from the 
BH (e.g. Cle net et al.l [2005'l or possibly stellar winds as 
discussed above. 

Based on the limited examples of M32 and 
NGC 404 (jSeth et al J 120101) . it appears NIR spectra at 
high spatial resolution may be an effective way to search 
for accretion in the lowest luminosity AGN. In both these 
cases, the NIR emission is the most luminous activity 
tracer, more luminous than the detected X-ray emission. 
Furthermore, the low luminosity of the X-ray emission 
(< io 38 ergs -1 ) ca n be confused with emission from 
X-ray binaries (e.g. iDesroches &: Hoi l2009t iGallo et al.l 
|2010|) and thus the NIR emission may represent the best 
indicator of BH accretion if it can be spectrally disentan- 
gled from the stellar light. One very interesting applica- 
tion of this technique would be to search for hot dust 
emission from accreting BHs in Gl and other globular 
clust e rs with putative c entral BHs (e.g. lUlvestad et al.l 
120071 INovola et alJl2010h . 

4. A KINEMATIC SIGNATURE OF PAST GAS ACCRETION 

In this section I discuss kinematic measurements from 
the CO bandhead of the NIFS data. The higher-order 
moments of the LOSVD show a very clear signature sim- 
ilar to that typically seen on larger scales in kinematic 
measurements of the centers of rotating, disky elliptical 



galaxies (e.g. iBender et a,l.lll994t IKrainovic et al.ll2008f ). 
The kinematic signature seen in these galaxies has been 
recently been simulated as resulting from gas accretion at 
the centers of g alaxies during the end phase s of gas-rich 
major mergers (Ho ffman et all [200l |20"Ioft . I analyze 
these kinematic results in light of these observations and 
simulations. 

Maps of the first four moments of the LOSVD derived 
from our NIFS data are shown in Fig. [2] The veloc- 
ity map shows strong rotation at a position angle of - 
25° (E of N) with an amplitude of ^55 kms -1 beyond 
a radius of 0'/3 (derived using the kinemetry software 
of IKrainovic et all 120061 ) . The dispersion increases to a 
maximum of 120 kms -1 in the central 0'.'5 due to the in- 
fluence of the BH and unresolved rotation; outside of this 
radius, it drops to values of ~60 kms -1 . The hi distri- 
bution is clearly anti-correlated with the velocity, while 
a more complex but axially symmetric pattern is visible 
in the M map with negative values along the minor axis 
and strong positive values along the major axis. The ef- 
fect of individual stars is clearly seen in the outer parts 
of the kinematic maps; these result in lower dispersions, 
extreme h3 values and strongly positive fe4 values. The 
prope rties of these stars are discussed in iDavidge et al.l 

mm . 

The top panels of Figure [3] shows the hZ-v/a and 
h4-v/a distributions of the data. The /i3 distribu- 
tion is strongly and tightly anti-correlated with v/a. 
The highest hA values are found in the most strongly 
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Fig. 3. — Distribution of the higher order LOSVD moments (/i3 and hi) as a function of the rotation velocity divided by the dispersion 
V/a at r < 1.5". The top panels display t he Ge mini/NIFS kinematics of M32, while the bottom panels display data from gas-rich mergers 
at similar orientations by Hoffman ct al. (2009). In the top panels, each spaxel in our IFU observations is plotted as a single data point 
(with a light gray error bar), while the simulation data has been analyzed in a similar way to mimic IFU data. 



rotating regions of the M32 nucleus. Both the /i3 
and M trends are remarkably similar to those seen on 
large r kiloparsec scales in observations (| Bender et al.l 
119941: iKrainovic et al.l [2008) and in recent simulations 
(Hoff man et al.ll2009[ ). The h3 anti-correlation with v/a 
is naturally explained by the presence of some stars cen- 
tered at the systemic velocity, while the positive M val- 
ues can also result from the superposition of a cold disk 
and hot bulge or be a signature of violent relaxation. In 
combination, these trends suggest the presence of two 
coincident components in the M32 nucleus, a dominant 
rotating disk and a thermalized, relaxed bulge structure. 

To facilitate visual comparison between our observa- 
tions and the simulations, in the bottom panels of Fig - 
ure [3] I have plotted models from lHoffman et al~l {2009). 
The plotted models are the results of 8 gas-rich (40% 
gas), equal mass mergers. The kinematic data was de- 
rived by viewing the central few kiloparsecs of these sim- 
ulations from a variety of lines-of-sight with inclinations 
consistent with th e 70±5° inclination fo r the central disk 
of M32 found bv iVerolme et al.1 (|2002l ). The slope and 
extent of the h3>-v/a distribution is almost identical to 
that seen in the simulations. For /i4, both the simula- 
tions and data have generally positive hi with a strong 
similarity in the shapes of the hA-v/a distribution. The 
very high h4 (> 0.15) values seen in the M32 data are 
the result of individual stars and clumps super-imposed 
on the velocity field. 



The kinematic data provide evidence for a two com- 
ponent system including a disk created by gas accretion 
into the nucleus. However, the details of the gas accretion 
leading to the formation of the central disk may be some- 
what different than in the simulations due to the large 
differences in scal e between them . Emb edded disk struc- 
tures occur in the I Hoffman et all (|2009f) simulations only 
in gas-rich major mergers and form on short timescales. 
However on the smaller scales of the M32 nucleus, it is 
possible that lower gas fraction mergers, gas-rich minor 
mergers, or even secular processes are still capable of cre- 
ating the disky kinematics seen in M32 on parsec scales. 
For instance, in iSeth et al.l (|2010fl . we present evidence 
that a minor merger ~1 Gyr ago in NGC 404 appears to 
be responsible for a burst of star formation that formed 
~5 million solar masses of stars in the central 10 pc of 
that SO galaxy. 

Another possibility for formation of nuclear disk is the 
gradual accretion of material from stellar winds at the 
center of the galaxy as described by |Baileyj (1980). This 
would nat urally lead to a gradua l trend in stellar age 
(jWorthevi [2001 iRose et all (20051) and the lack of an y 
break in the rotation curve ([Dressier fc Richstonel [l988). 
Furthermore, this scenario may nat urally account fo r the 
abundance gradie nt in the nucleus; iWorthevI ([2004D and 
IRose et all pOOl find the [Mg/Fe] ratio is significantly 
subsolar in the nucleus and increases with increasing ra- 
dius. If the stars in the nucleus were made from the 
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stellar winds of previous generations of stars, then this 
might result in a drop in the [Mg/F e] (e.g. as in the Milk y 
Way globular cluster NGC 2808; iCarretta et al.l [2009h . 
There is sufficient mass in the bulge of M32 to supply 
the necessary gas to the nucleus. However, whether the 
gas would actually end up as a rotating structure at the 
center would require detailed hydrodynamic simulations. 
This scenario could also be tested using stellar popula- 
tion modeling of resolved stars and spectra in and around 
the nucleus if the models are accurate enough to distin- 
guis h between populatio ns of intermediate (2-10 Gyr) age 
(e.g. iMarigo et aLll2008| ) . 

4.1. Asymmetries in the M32 Nucleus 

Simulations by iHopkins fc Quata"ert (2010a. b) follow 
gas accretion down to parsec scales in galaxy mergers and 
find that gas accretion in the presence of BHs cause the 
formation of lopsided disks. These disks corr espond quite 
well to the eccentric disk seen in M31 (e .g. lLauer et al.l 
119931) . and IHopkins fc Quataertl (|2010bf ) suggest that 
this disk formed at the same time the M31 BH was grow- 
ing. Due to the parsec scale of the expected asymmetries, 
M32 is one of the few other galaxies close enough for the 
asymmetries to be seen. If the M32 nucleus had a ma- 
jor mfir^e^JojnnaHpji_^enajipsimilar to that considered 
by IHopkins fc Quataertl (|2010b|) . a lopsided nuclear disk 
would be expected. While M32 clearly does not have the 
prominent asymmetric structure similar to that seen at 
the center of M31, I analyzed HST images to see if there 
is any evidence for a lower level asymmetry. 

Archival HST WFPC2/PC images were obtained in 
four bands (F336W, F555W, F675W, and F814W) from 
the CADC WFPC2 association^. I then fit axially sym- 
metric model profiles to the data, subtracted these off 
and examined the residuals. Fits to the images in each 
filter were performed on the central 20" of each image 
using elliptical double Sersic profiles. The outer Sersic 
component was constrained to have the same indices and 
effective r adii as fi t s to t he 1- D M32 surface-brightness 
profile bv lGrahaml ([2001 and lGraham fc Spitleri (2009) 
that include measurements from the same F814W HST 
data we use here. All fits were c onvolved wi th an ap- 
propriate TinyTim WFPC2 PSF (|Kristlll995| ). The re- 
sulting fits accurately described the light profile at all 
radii to within 3%. The inner Sersic profile fits had ef- 
fective radii of 1.5 - 1.6", consistent with that found by 
IGraham fc Spitlerl ff2009h . minor-to-major axis ratios of 
0.73 and position angles of —20°. After subtraction of 
these profiles, some asymmetric residuals were visible 
(see Fig. |4j left panel) . To quantify these, I compared 
the total residual flux on one side of the minor axis to 
the other in elliptical annuli (Fig. @] right panel). At 
radii between 0'.'45 and 1"8 (~2-7 pc) there is an excess 
of 1-2% of the total flux in the southern half relative to 
the northern half. This is seen most clearly in F555W 
and F675W, and is also present but less clear in the other 
two bands due to lower S/N (at F336W) and increased 
surface brightness fluctuations (in F814W). 

This small photometric asymmetry is much less lop- 
sided than the central part of M31 (where the disk ec- 
centricity ~40%). The asymmetry corresponds to a total 
mass of ~4x 10 5 M derived from the F814W band using 



the mass-to- light ratio from iVerolme et a l. (2002). It is 
unclear if such small a symme tries are consistent with the 
IHopkins &; Quataerti (|2010b[ ) model. The physical situ- 
ation in M32 is quite different from M31, with a much 
smaller BH mass and a higher stellar density. Simula- 
tions that are more appropriate to M32 may serve to 
both (1) help understand in detail how the M32 nucleus 
formed, and (2) test the feasibility of the lopsided disk 
model for driving BH growth in lower mass systems. 

5. CONCLUSIONS 

In this paper, I have presented results on the current 
and past gas accretion in the M32 nucleus using adap- 
tive optics assisted Gemini/NIFS data. The data reveal 
an unresolved emission component at the center of M32 
that results in the weakening of CO lines and redden- 
ing of the continuum slope in the nucleus. This emis- 
sion appears to originate from hot dust emission close 
to the BH as seen in many higher luminosity AGN. The 
hot dust emission is unresolved with a size of <0.9 pc 
and has a luminosity of ^2xl0 38 ergs -1 . This luminos- 
ity is more than two orders of magnitude more luminous 
than the X-ray emission detected from the central BH by 
IHo et al] ()2003[ ) , the only other bandpass in which emis- 
sion due to accretion has been detected in this galaxy. 
Combined with similar emission in the nearby low lumi- 
nosity LINER NGC 404, this suggests that high resolu- 
tion NIR spectra may be an efficient way of locating low 
level activity from BHs in nearby galaxies and globular 
clusters. 

I also present the highest quality kinematics of the M32 
nucleus currently available, potentially useful for any fu- 
ture dynamical studies of this benchmark system. These 
measurements cover the central 3" x 3" with a resolution 
of 0"25. The very high signal-to-noise of these measure- 
ments allows accurate determination of the higher order 
moments which reveal a remarkable kinematic signature 
similar to those seen on larger scales in elliptical galax- 
ies and merger simulations. The signature manifests as 
an anti-correlation of ft.3 with v/a and strongly positive 
hA especially in areas of strong rotation. These kinemat- 
ics suggest the presence of a dominant disk component 
embedded in a pressure supported system. I propose 
a possible formation scenario in which the M32 nucleus 
forms from stellar winds of the M32 bulge. This sce- 
nario qualitatively explains the kinematic, stellar pop- 
ulation and abundance properties of the nucleus. I 
also briefly consider whether M32 contains any asym- 
metries in its nucleus , as mig ht be expected based on the 
Hopk ins fc Quataerti (|2010bD simulations of major merg- 
ers. A small ~2% photometric asymmetry is found at 
radii between 2 and 7 pc. 
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Fig. 4. — Figures showing a small asymmetry in the M32 nucleus. Left - shows the a model-subtracted F555W image oriented with N 
up and E to the left. Black solid lines show the annulus with major-axes of 0.45" to 1.8", while the dashed line separates the northern 
and southern halves of this annulus. Right - comparison of the residuals of the model-subtracted images. The Y-axis plots the difference 
in flux between the southern and northern annuli in the model-subtracted images normalized by the original image. A 1-2% excess in the 
southern annulus is seen in all bands from 0.45— 1.8" (2—7 pc). 
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